Central-moment description of polarization for quantum states of light 
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We present a moment expansion method for the systematic characterization of the polarization properties of 
quantum states of light. Specifically, we link the method to the measurements of the Stokes operator in different 
directions on the Poincare sphere, and provide a method of polarization tomography without resorting to full 
state tomography. We apply these ideas to the experimental first- and second-order polarization characterization 
of some two-photon quantum states. In addition, we show that there are classes of states whose polarization 
characteristics are dominated not by their first-order moments (i.e., the Stokes vector) but by higher-order pola- 
rization moments. 

PACS numbers: 42.25.Ja,42.50.Ar,03.65.-w 



I. INTRODUCTION 



under any polarization rotation. 



A fundamental property of light is its vector nature. Far 
from a source, freely propagating light can be approximated 
by a plane wave, with the electric field directed in the plane 
perpendicular to the direction of propagation. Already the 
early pioneers of optics realized that a convenient way of char- 
acterizing light is to describe the figure the tip of the electric- 
field vector traces out in this plane. Stokes established an op- 
erational procedure to characterize not only the polarization 
properties of light, but also to what extent a field is polar- 
ized IJJ. The method is still the standard way of assessing 
polarization, although several generalizations, such as polari- 
zation of non-plane |2-5| and multi-mode |6-8| fields, have 
been developed. A limitation of Stokes' approach is that it 
only considers the average intensities (or photon numbers) and 
hence only assesses the first-order polarization moments. 

As polarization is a relatively robust degree of freedom, 
that, moreover, can almost losslessly, cheaply, and easily be 
transformed, it is very often used for coding and manipu- 
lating quantum information. Examples of experiments rely- 
ing on polarization include quantum key distribution ll9l [T0l . 
quantum dense coding iTT\, quantum teleportation lfT2l . quan- 
tum tomography llT3l . rotationally invariant states lfT4l . phase 
super-resolution ifTSi , and weak measurements fT6l . How- 
ever, many of these experiments use correlation measure- 
ments, effectively using second, or higher, polarization mo- 
ments. Such correlation measurements can give surprising 
results. For example, states that appear unpolarized (that 
is, with vanishing Stokes parameters), can show unit visi- 
bility polarisation-correlations when rotated on the Poincare 
sphere [17J. Such states have been said to have "hidden pola- 
rization" lfT8l[T9l . As we shall discuss below, there are actu- 
ally large classes of such states, and they can be classified by 
the number of lowest-order moments that are invariant under 
polarization transformations. We shall refer to such states as 
rth-order unpolarized if the first rth moments are all invariant 



As hinted by this discussion, the full description of polari- 
zation can be sorted into moment orders, and simultaneously 
(but perhaps less obviously) into excitation manifolds. A con- 
venient and experimentally palatable way to do this is by the 
use of central moments. 

For the three lowest orders, the central moments coincide 
with the cumulants introduced by Thiele |20|. Each succes- 
sive cumulant provides information of statistics not already 
contained in the lower-order cumulants. They have some 
advantages over a moment description when making affine 
transformations, and they also provide a simple method of 
quantifying the difference between a statistical distribution 
and its simplest Gaussian approximation 11211 . (For Gaus- 
sian distributions all cumulants of order > 3 vanish.) Kubo 
promoted their use in quantum mechanics and thermodynam- 
ics 1221 . but in polarization optics they have been used rather 
sparsely 123^271 . 

Below, we first recall some definitions and notation in 
Sec.|ll] In Sees. UlllandllVlwe examine how first- and second- 
order polarization properties can be described in terms of ex- 
pectation values and central moments, respectively. In the fol- 
lowing two sections, [V]and [Vl] we discuss how the formalism 
can be extended to orders higher than the second. In Sec. VII 



we subsequently discuss the connection between excitation 
manifolds and polarization data and show that polarization to- 
mography in general requires far less data than full state to- 
mography. We then apply the formalism, both theoretically 



and experimentally, to certain polarization states in Sec. VIII 



In particular, we show that there are many states whose pola- 
rization characteristics are dominated not by their first-order 
moments (i.e., their Stokes vector), but by higher-order pola- 
rization moments. For example, for three-photon states there 
exist six different classes of states with different polarization 
characteristics. Finally we draw some conclusions from the 



analysis in Sec. IX 
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II. STOKES OPERATORS AND THE STOKES VECTOR 

We will build on the classical theory of polarization based 
on the Stokes parameters. For quantized fields, the Stokes 
operators [28 1 take the role of the Stokes parameters. They 



^2 — i{d}]dy — flj^ay) , ^3 = ci^^d}{ — ayCiv , 



(1) 



where uh and ay are the annihilation operators of the two or- 
thogonal modes, in the following taken to be linearly hori- 
zontally and vertically oscillating electric fields, respectively. 
The annihilation operators obey the bosonic commutation re- 
lations 



= 5 



j,ke{H,v}. 



(2) 



The average values of the Stokes operators correspond to 
the Stokes parameters ({So), (S)), where the Stokes vector op- 
erator S is S = (Si,S2,Sj,). In terms of the Poincare sphere, 
the definitions mean that ^2 is the eigenoperator for a cir- 
cularly polarized field, and thus that the operator is parallel 
to the axis through the south (left-handed circular) and north 
pole (right-handed circular) of the sphere. Si and ^3 are the 
eigenoperators for diagonal and anti-diagonal, and horizontal 
and vertical, linear polarization, respectively. These operators 
lie in the equatorial plane of the Poincare sphere. The direc- 
tions of §1, §2, and ^3 form a right-handed orthogonal vector 
set in the Poincare space. 

The Stokes operators satisfy the commutation relations of 
the su(2) algebra; 



[Sj,Sk] — '2£jkeSe , 



where Ejiic is the Levi-Civita tensor The non-commuting 
character of these operators leads to the uncertainty relation 



2(5o)<(S2)-(S)2<(5o>((5o) 



(4) 



In spherical coordinates we can use the polar and azimuthal 
coordinates 6 and (p to parameterize the unit vector as n = 
(sin cos 0, sin sin 0, cos 0). (Note, however, that is the 
angle to ^3, and that §1 and ^3 lie in the equatorial plane of 
the Poincare sphere, as explained above.) We can now express 
the Stokes operator in any direction n as 



= S • n = niSi + 11282 + n^S3. 



(5) 



In addition to the commutation relation Q, one also has the 
relation 



[So,Sj]=0, ./e {1,2,3}. 



(6) 



This indicates that there exist simultaneous eigenstates of Sq 
(giving the total photon number) and any other Stokes opera- 
tor So in principle, a measurement of §„, if repeated on many 
members of an identically prepared ensemble, also allows the 
photon number statistics to be determined. In fact, a common 



way to measure any Stokes operator §„ ~ UaS^U^, where Ua 
is a (unitary) linear polarization transformation rotating the 
axis 3 to align with n, is to first "rotate" the state according 
to p UnPU„, and then measure ^3. That is, after the rota- 
tion of the state, one separates the H and V modes by polari- 
zation optics and then counts the number of photons in each 
mode. The photo-count difference then gives the measured 
5n eigenvalue, while the sum gives the Sq eigenvalue. This 
suggests that in a full description of quantum polarization, the 
excitation manifolds should be treated separately 129^-^311. In 
consequence, coherences between different manifolds do not 
carry polarization information. Below we shall use the total 
photon number as an index of the excitation manifold. As 
it will simplify the subsequent discussion, we introduce the 
normalized A^-photon density matrix defined as 

Pmn.N = —{m,N -m\p\n,N -n), OT,n £ {0, . . . ,A^}, (7) 
Pn 

where = Y!^^^Q{m,N ~m\p\m,N — m) . With this definition, 
we have 



Tr(p^S). 



(8) 



In reality, it may be experimentally difficult to divide the po- 
larization measures into excitation manifolds, except for few- 
photon states. To this end, we shall also define the excitation 
averaged Stokes vector 



(S)= Lp^(S) 



(9) 



^=1 



All other measures of polarization, defined below, can be av- 



(3) eraged over the manifolds in the same manner In Sees. VIII C 



and VIII D this is done when discussing two-mode coherent 
states and two-mode thermal states, states that when they are 
used often contains a large, but indeterminate, number of pho- 
tons. However, for few photon states it is possible to divide 
the results according to photon number through coincidence 
measurements, and we believe that it soon will become com- 
mon to use detectors with photon number resolving capability 
for such states. 

The idea that we will develop below is that the rth-order 
polarization in the A^th excitation manifold is characterized by 
a data set that can predict {S^)n for any direction of the unit 
vector n on the Poincare sphere. 



ni. FIRST-ORDER POLARIZATION MOMENTS 

Since the classical description of polarization is based on 
the first-order moments, the quantum description is the direct 
translation of the classical description. That is, the Stokes vec- 
tors (S) N defined in (jsj gives the complete first moment pola- 
rization information. It follows from the expectation value of 
both sides of Eq. (jsjl with regards to the state pj^, that (S) n is 
sufficient to predict Tr(j3Af5n) for any n. 
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IV. ASSESSING THE SECOND-ORDER POLARIZATION 
MOMENTS 

How should one then go about to characterize higher-order 
polarization properties? One way would be to assess all 
second-order moments, i.e., all polarization correlation values 

of the form r^^'^' (p) =Tr(pA,5j4),;,'te {1,2,3}. However, 
only when j = k these operator products are Hermitian, so the 
expectation values cannot be measured directly. Nonetheless, 
from a theoretical perspective such an approach is viable and 
equivalent to the description via polarization moments in dif- 
ferent directions. In |32| we have followed this path. A great 
simplification and reduction in data is to collect the polariza- 
tion correlation information into Hermitian moment compo- 
nents |32|. One of this method's advantage is its simple hier- 
archy over the moment orders. It is straightforward to under- 
stand how to systematically collect the needed, non-redundant 
information moment order by moment order. Experimentally, 
it is equivalent to the proposed method in that one measures 
successively higher moments of the Stokes operator for se- 
lected directions on the Poincare sphere and then solves an 
ensuing equation system. A drawback is that it is not so easy 
to see the relative importance of the moment orders, as lower 
order moments contribute to higher order moments. 

Another method is via the two-mode coherence matri- 
ces lfT9l where the rth order coherence matrix coefficient jk 
is defined by {{ajj)-' {ayY^ ^ a'^a'^'') . This method is informa- 
tionally equivalent to the method we shall develop. Among 
the advantages with this method is that all moment coeffi- 
cients are expectation values of normally ordered annihilation 
and creation operators, making calculations for coherent states 
particularly easy. Another advantage is that for -photon 
states, all coherence matrices of order r > N vanish. How- 
ever, the method has only an indirect connection to the Stokes 
operators, and the ensuing matrices give little direct "feeling" 
for the polarization properties of the state, although they con- 
tain all the needed data. Experimentally, the off-diagonal co- 
efficients j k of the coherence matrices are not straightfor- 
ward to measure. In ||33l a method using phase plates and 
projection onto an rth order "intensity" of one of the modes, 
namely ((a^)''a^), is proposed (note that this is not equal to 
a measurement of ((aj^o//)'')). Choosing properly {r + 1)^ 
different settings of the phase plates and solving the ensuing 
set of linear equations the coherence matrix of order r can be 
obtained. For a state containing up to, and including, pho- 
tons, A^(2A^^ +9A^+ 13) /6 measurements are thus required us- 
ing the measurement scheme proposed in 1331 . roughly twice 
more than for our scheme, see the end of Sec. |VI| and Sec. |VII| 
below. The proposal in |33| discuss only characterization of 
A^-photon states, and not how to assess the polarization or 
higher order coherence-properties of states with an indeter- 
minate number of photons. 

To see how the polarization central moments appear quite 
naturally in a polarization description, we expand each opera- 
tor in a state-dependent mean and a fluctuation part, v.i.z. 



In the following, to simplify the notation, we shall write 

Tr(pw50 = {S'n)N and Ti[pnKAp)] = {K)n- This aUows 
us to write, for r = 2 



n]{{Si)l + {A\)N)+cycl. + cycl. 

+nin2{2{Sy)N{S2)N + (AiA2>7V + (A2Ai)Ar) 

+cycl. 

+cycl., (11) 



where cycl. denote a cyclic permutation of the indices. We 
see that apart from {Sj)^, j E {1,2,3}, the expectation values 
of the six Hermitian fluctuation "operators" in ( 11 1 are what is 
needed to know (5^) n in any direction. These expectation val- 
ues are the second-order central-moments, (coinciding with 
the second-order cumulant) defined as 



(A/A<.)Af — {SjSi.)N - {Sj)N{Sk)N- 



(12) 



As can be seen from (H) it is convenient and natural to col- 
lect the mixed-product {j ^ k) central moments into Hermi- 
tian terms, e.g., (AyA^. + A<.Ay)Af. These terms can be mea- 
sured, and we see that in addition to the Stokes parameters, 
we need six more numbers to fully characterize the second- 
order polarization-properties. The first three can be obtained 
from measuring the statistics of the Stokes vector S yield- 
ing the first-order moments {{Si,S2:St,))n and the variances 
{A^)n7 j G {lj2,3}. The additional three numbers can be ob- 
tained from measuring the statistics of Sn along the "diago- 
nal" directions (l,l,0)/%/2, (l,0,l)/%/2, (0, 1, l)/%/2 in the 
S1S2, SiSi, and ^2^3 planes, respectively, corresponding to 
the angles (0,0) of (7r/2, 7r/4), (7r/4,0), and (7r/4,7r/2) on 
the Poincare sphere, and then using ( [TT| . 

As a minor digression, these second-order central-moment 
terms are directly connected to the Hermitian polarization co- 
variance matrix Tn with matrix coefficients 



'^jk.N -{AjAk + ^k^j)N, 



(13) 



where j,k E {1,2,3} ||34]| . Each such matrix has six indepen- 
dent coefficients as ^jk.^^ = ^kj,N by construction. It is clear 
from Eq. 
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An,iv(p) = 5n - Tv{pNSn) 



(10) 



that this covariance matrix contains the infor- 
mation we need, in addition to the expectation value of the 
Stokes vector, to be able to predict the value of {S-^)n in any 
direction. We also have {A^)n = n • T^; • n', where t denotes 
the transpose. 

Every covariance matrix Tn can be made diagonal by an 
orthogonal matrix R. In this rotated, orthogonal coordinate 
system, where Sej point in the direction of eigenvector e/, j E 

{1,2,3} of Fa?, one finds the extreme values of {A'^)n- In this 
coordinate system Eq. ( [TT| i simplifies to 

(A^> N = ^i (sin 0' cos (j)')^ + A2(sin 0' sin + A3 cos^ 0', 

(14) 

where Xj is the jth eigenvalue of Tn, 9' is the angle between 
n and 63, and (p' is the azimuthal angle in the 61-62 plane. This 
equation may look like the equation of an ellipsoid, but it is 
not, as this is the magnitude of the variance of (A^) in the 
direction n on the Poincare sphere. 
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In order to measure Tn, one makes the same measurements 
as were discussed above. The matrix Tn can subsequently 
be deduced by solving Eq. (jTT| for {Sj)^ and {A^j)n given 

the measured values of {S^)n along the six directions. For 
better "immunity" to systematic errors, one could make mea- 
surements along additional directions and subsequently make 
a best fit of the ensuing overcomplete system of equations. 



V. THIRD-ORDER POLARIZATION 

Moving on to third-order moments, things get a bit more 
involved. Still, our underlying idea is that if one has all the 
central moments up to order three, then one can predict {S^) 
for any direction. 

We therefore first express the expectation value (5n)^ in 
terms of (S)n- 



{Sn)N 



+3{nln2{Si)jj{S2)N + nini {Si)jj{Si)N 
+n2n\ {S2)1,{S\)n + n2n^{S2)l!{Si)N 
+n\n\ {Si)]^{S\)N + n],n2{S-i)l^{S2)N) 
+6nin2n3{Si)N{S2)N{S3)N- (15) 



In a similar manner we can express the third-order raw mo- 
ment of Sn as 



{sI)n 



2\ 

i)n 



+cycl. +cycl. 

+«]n2(3(5i)^(AiA2 + A2Ai)a, 
+3(52)7v(A?)iv + (A?A2+A2A 
+ (AiA2Ai)A,)+cycl.+cycl. 
+«?«3(3(5i)iv(AiA3+A3Ai)A, 
+3(53);v(A?)iv + (AlA3+A3A2);v 
+ (AiA3Ai)A,)+cykl.+cykl. 

+nin2n3(3(5i)iv(A2A3 +A3A2)iv 

+ (AiA2A3+AiA3A2)Af 

+cycI.+cycl. + (5„>^). 



(16) 



Finally we can express the third-order central-moments as 

{AjAkAi)^ = {SjSkSe)N - {Sj)N{SkS(')N 

-{Sk)N{SjSi)N ~ {Se)N{SjSk)N 

+2{Sj)N{Sk)N{Si)N- (17) 

Hence, for the first to third order, the central moments coin- 
cide with the cumulants. 

One sees that in ([T6|, if the ten Hermitian, third order, cen- 
tral moment terms, each associated with a different geometric 
term njiTkn^^j^k^ G {1,2,3}, y' + A: < 3 are determined, 
in addition to the first and second-order properties, then the 
third-order polarization-properties are also determined for any 
direction. Hence, what one needs to measure are the sums 



of all fluctuation terms having j ones, k twos, and 3 — j ~k 
threes, where j + k < 3, or more generally, for polarization 
order r, into sums having r — j — k threes, where j + k < r. 

Measuring the third-order fluctuations along, e.g., 
the (0,0) directions (0,0), (7r/2,0), (;r/2,7r/2), 
(7r/2,(^)i), (7r/2,-(/)i), (jt/l - 0i,O), (7r/2 + 0i,O), 
{n/2 - (j)un/2), {n/2 + ^i,k/2), and (7i:/2 - 0i,7r/4), 
where 0i = arccos \/2/3, one gets a system of ten linearly 
independent equations that allows one to determine the terms 
(Aj)ftr + 3(5i)Af(A[)A? etc. Using the knowledge about the 
lower-order polarization terms, one can subsequently estimate 
the third-order terms, in this case {A\)n- We note that the 
three first measurement directions are simply along the Si, 
§2 and ^3 axes, so in fact, only measurement along seven 
extra directions are needed, in addition to the measurements 
along six directions needed to determine {S)n and {S-^)n- 
Alternatively, if one wants to minimize the number of mea- 
surement directions, one can use the statistics collected when 
measuring along the six directions that determine the first and 
second-order polarization moments, and supplement them 
with measurements along the four new directions (;r/6, 7r/6), 
(;r/6,;r/3), (;r/3,;r/6), (;r/3,;r/3). 

For third-order polarization the first thing to be considered 
is that the fluctuations of involves not only third powers 
of Aj, but also terms like n\{Si) {Af)^ and n\n2{Si)N{AiA2 + 
A2Ai)Ar. That is, the second and the third-order fluctuations 
become "intermixed" in this polarization order unless the state 
has vanishing Stokes parameters. This is in contrast to the 
(simpler) second order A consequence of this is that if the 
state's first-order polarization is much larger than the square 
root of its variance, then all higher-order fluctuations will, in 
general, be dominated by the beating terms between the mean 
polarization vector and the second-order fluctuations. Hence, 
for most "reasonably excited" and "somewhat first-order po- 
larized" states one needs not go beyond the second-order mo- 
ments to characterize the polarization fluctuations of all or- 
ders to a very good precision. However, for states having a 
small or vanishing first order polarization, and for, e.g., the 
eigenstates to the Stokes operators in the direction of Fat's 
eigenvector directions on the Poincare sphere, the polarization 
structures of orders higher than two will be of interest. 



The expansions (111 and (16 1 also indicate an experimental 
advantage in describing the polarization in terms of increasing 
orders of its central moments. For each order it becomes quite 
clear to which accuracy one needs to measure the moments 
to obtain information not already contained in lower moments 
and similarly, to what extent the higher-order central moments 
contribute to the raw moments. This information is of course 
implicit in "equivalent" descriptions such as generalized co- 
herence matrices 1 19] or polarization tensors l,32J . but it is not 
explicitly displayed. 



VI. FOURTH- AND HIGHER-ORDER POLARIZATION 

From the preceding sections it is rather clear how one could 
continue through the higher orders. In order to know {S'^^n in 
any direction, the full set of Hermitian central-moment terms 
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for all orders < r is needed. 



Explicitly, the fourth order the central moment is 



(AyAiA^A,„)Af = {SjSkSeS,„)N — {Sj)N{SkS(Si„)N — {Sk}N{SjS(S„j)N— {S()N{SjSkS,„}N 

^{Sm)N{SjSkSi)N+ {Sj)N{Sk)N{SeS,n)N+ {Sj)N{Se)N{SkSm)N + {Sj)N{Sm)N{SkSe)N 
+ {Sk)N{Se)N{SjSm)N + {Sk)N{Sm)N{SjSe)N + {Se)N{S,„)N{SjSk)N 
— 3{Sj)N{Sk}N{S(:}N{Sm)N- 

I 



(18) 



In contrast to the three lower orders, this result is not identical 
to the fourth-order cumulant. Higher-order central moments, 
that we will not write out explicitly, do not coincide with the 
cumulants either. 

In analogy with the second and third order, we need not 
determine each term of the form ( 18 1, but only the Hermitian 



sum of moments associated to a certain geometrical pre-factor 
The number of such central-moment sum-terms specific to the 
order r is (r + 1 ) (r + 2)/2 and the complete set of such terms 
up to, and including, order r is r(r^ + 6r + ll)/6. To obtain the 
terms, one would have to measure the polarization statistics 
along such a number of carefully selected directions, yield- 
ing a complete, "maximally" linearly independent set of equa- 
tions that could be solved numerically. To obtain better ac- 
curacy one could "oversample" the polarization statistics over 
the Poincare sphere and use maximum likelihood or entropy 
methods to make a better estimate. However, as the states that 
have their polarization characteristics mainly determined by 
the rth-order moment will be rather elaborated as r increases, 
the interest in the polarization central moment terms will be 
limited to r < 4, or so. 



VII. POLARIZATION PROPERTIES AND EXCITATION 
MANIFOLDS 

Using the bosonic commutation relation it is possible to 
rewrite any rth-order product of Stokes operators to a sum of 
normally ordered creation and annihilation operators of max- 
imum order r in the annihilation orders II33L To exemplify, 
one can write 



Si 



-\-aJjCiff - 



' ayUy . 



(19) 



As all Stokes operators are composed of terms with one cre- 
ation and one annihilation operator, this implies that all pola- 
rization properties of a state with no excitation above the A^- 
photon manifold are determined by the polarization moments 
up the r = Nth order. All moments of order higher than 
can have only those normal ordered terms less or equal to the 
A^th order different from zero, and those terms will always be 
contained in the moments up to, and including, the A^th order 
Below we shall see a specific example of this, namely that 
for a three-photon state, it is sufficient to require that (5„ ) is 



isotropic for m = 1,2,3 in order for the state to be unpolarized 
to all orders. Note, however, that should the higher-order cen- 
tral moments be zero, this does not indicate that the state lacks 
higher-order polarization-structure. Instead, the implication 
is that this structure can be derived from the "beating" terms 
from lower-order polarizationmoments, as already hinted in 
Sec.|V] 

Another consequence of the fact that states with no exci- 
tation above the A^-photon manifold has its polarization fully 
characterized by its A^ lowest-order moments is that polari- 
zation tomography of such a state is requiring considerably 
less resources than a full state tomography. For a full state 
tomography involving the (A^+ l)(A^ + 2)/2 basis states (e.g., 
for A^ = 1 the states |0,0), |0, 1), and |1,0) can be chosen) the 
density matrix is characterized hy N{N^ +6N^ + 13A^+12)/4 
independent real numbers. This can be compared to the 
A^(A^^ + 6A^ + 1 1)/6 numbers needed for the polarization to- 
mography of such a state. Raymer et al. has used the term 
"polarization sector" of the density matrix for the subset of 
information needed to characterize only a state's polariza- 
tion m. 

This said, an A^ -photon state is fully described by A^(A^ + 2) 
real numbers while the polarization central moments up to, 
and including the r = A^:th order require 
numbers. Very recently we found a method to reduce the 
number of measurements to the minimum N{N + 2), but the 
analysis and description of this scheme will be published else- 
where li32J. 



VIII. APPLICATION TO DIFFERENT POLARIZATION 
STATES 

We shall now apply the characterization developed above to 
a few examples and also compare the theory with experiments 
in the two-photon excitation manifold. We remind the reader 
that we use the ^3 eigenstates as our basis states. The experi- 
mental setup is discussed, and measurement data are given, in 
the Appendix. 



A. SU(2) coherent states 

Through an appropriate polarization transformation of the 
state \N)h'^\0)v = \N,0) any A^-photon, SU(2) coherent 



-r~----.,_ ; / i / FIG. 2: (A^) in the plane defined by (A^) = for the state |2,0). 

^^2^^ ""p-:'' ^ '^j'-i'' The theoretically expected result is drawn solid (blue), the experi- 

mentally obtained result is drawn dashed (red), and the experimental 
result, solidly rotated so that its eigenvectors coincide with the in- 
FIG. 1: The second order central moment (A„) for the state |2,0). j^^^^^ eigenvectors, is drawn dotted (black). 
Theoretical plot in (a) and the experimental results in (b). 



state can be obtained. Since a polarization transformation 
is equivalent to a rotation of the Poincare sphere, it thus 
suffices to study the state |A^,0). Quite clearly, all its mo- 
ments are zero except in excitation manifold A^, and there- 
fore we will suppress this index. The state has (So) = N, the 
Stokes vector is (0,0, A?), (A'^') = V m, (Af) = (A^') = 
for odd m, (A^) = (A^) = A^, and (A1A2 + A2A1) = (A1A3 + 
A3A1) = (A2A3 + A3A2) = 0. Its second-order, polariza- 
tion central moments are hence reduced to a toroidal struc- 
ture with radius A^, with its "hole" in the ^3 direction on 
the Poincare sphere. Its third-order central moments have 
the non-vanishing terms (AyA3 + AyA3Ay + A3A^) — —2N, 
j E {1,2}. The non-vanishing terms of fourth order are 
(Aj) = 3N^ - 2N, {AjAl + A^A] + A;A3 A^'A3 + A3 AyA3 A,- + 



A.A^A, 



- A3AjA3 

{2 



4N, and (Af A| + A^Aj 



{2X2 _ 



A2A1A2A1 +AiA^Ai + A2Af A2) = 6A^- 



A1A2A1A2 + 
4A?,;e {1,2}. This 



is a minimum-sum uncertainty- state [saturating the left in- 
equality in Eq. Q]. In Fig. [T]we plot the theoretically com- 
puted function (A^)2 for the state |2,0) to the left, and the 
experimentally obtained results on the right. The measured 
Stokes vector of this state is (—0.19,0.12,1.97). The main 
source of error in the estimation of the Stokes vector is neither 
fluctuations nor random errors, but the fact that the generated 
state and the measurement axes are slightly rotated relative 
to each other. This is better seen in Fig. [2] where we have 
plotted the theoretically expected result (solid) and the exper- 
imental figure (dashed) as deiived directly from the measured 
parameters in Table |ll| (no fitting done) in the §1 -Sj, plane. 
We have subsequently used the measurement data and com- 
puted the eigenvectors and eigenvalues for the matrix r2 de- 



fined in Eq. ( 13 1. The eigenvalues (that due to the Hermitic- 
ity of T2 are the plot's extrema which are found at orthogo- 
nal points on the Poincare sphere) are 2.08, 2.02, and 0.00, 
respectively, and, e.g., the eigenvector coiTesponding to the 
smallest eigenvalue is (—0.15,-0.02,0.99). This eigenvector 
is rotated about 8.1 degrees from the ^3 axis, mainly around 
the ^2 axis (roughly in agreement with the orientation of the 
state's Stokes vector). If we re-plot the figure, with no other 



"fitting" than a solid rotation of the experimental figure to 
make its eigenvectors coincide with the Poincare sphere co- 
ordinate axes, we obtain the dotted curve in Fig. |2] Note that 
the needed rotation is not exactly perpendicular to the drawn 
plane, so the rotated figure will also change shape slightly. 
The dotted figure, (and the values of the T2 eigenvalues) con- 
firm that the experimental errors due to fluctuations are below 
±4%. The main errors are systematic, due to imperfect pola- 
rization optics and beam splitters. The systematic errors will 
persist even if additional points on the Poincare sphere are 
measured, indicating that the proposed method is "efficient" 
from a data collecting point of view. By measuring the nine 
data in Table [III one ^'^^ hence estimate both the accuracy and 
the precision of the experiment. 



B. |A',A^) states 

This state has {§0) = 2A^, the Stokes vector is (0,0,0), 
and (A'3 ) = V m. The only non-vanishing second-order, 
central-moment terms are (A?) = 2A^(A^+ 1), j G {1,2}. The 
state has vanishing third-order central-moment in every di- 
rection, and its fourth-order, non-vanishing central-moment 
terms are (A;*) = 2A^(3A^^ +6N^+N- 2), (A^A^ + A^A^ + 

A;A3A^A3 + A3A;A3A^' + A^A^A^ + A3A2A3> ^ m{N + i), 

and (A^A^ + A^A^ + A1A2A1A2 + A2A1A2A1 + AiA^Ai + 
A2A2A2) = 4N{3N^ +6N^+N- 2), ./ G {1,2}. This is a 
pure, maximum-uncertainty state, saturating the right inequal- 
ity in Eq. (|4|i]. If one plots the experimentally obtained re- 
sults one obtains a figure very similar to Fig. [T] (b), except 
that (A^)2 for this state is approximately twice as large in all 
directions as for the state |2,0). The measured Stokes vec- 
tor of this state is (-0.01,-0.08,0.01). Extracting the mea- 
sured state's extremal values for (A^)2 from the experimental 
data (i.e. the eigenvalues of r2) one obtains 4.10, 3.98, and 
-0.03, indicating an experimental error due to fluctuations be- 
low ±3%. This data too, indicates a slight mismatch (about 
12 degrees) between the state's intended orientation on the 
Poincare sphere and its measured orientation. Note that the 
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state's orientation cannot be obtained from the Stokes vector, 
since nominally it vanishes. The the measured values are sim- 
ply the measurement errors. However, from the second-order 
moment data the state's orientation relative to the Poincare 
sphere axes can be well determined. 

C. Two-mode coherent states 

Any two-mode, coherent state \a' .a") can be converted 
into the state ||o:|,0), where |ap = |a'p + \a"\~, by a polari- 
zation transformation. Therefore it suffices to study the latter 
state, which can be written 

exp(-|a|V2)f ^|A^,0). (20) 

In each excitation manifold except the non-excited one, the 
state has the same central moments as a SU(2) coherent 
state. Summing over the manifolds, the coherent state has 
(5o) = |ap, the Stokes vector (0,0, |ap) and (A^) = |a|2 
for j G {1,2,3}. The off-diagonal coefficients of the covari- 
ance matrix T are zero, so the second-order central moment 
is isotropic with radius |ap. In the third order, the only 
non-vanishing central-moment terms are: (A3) = (AjAs + 
AiAsAi+AsAf) = (A^A3+A2A3A2+A3A^) = |a|2. The non- 
vanishing fourth-order central-moment terms are: 

(A))=3|ar + |a|2, 

for i G {1,2,3}, and 

{k]kl + kjhkAk + kjklkj + hk]kk 

+kkkjkkkj + klk^j)^6\a\'^ + 2\a\', 
for 7, ^ e { 1 , 2, 3} and j < k. 

D. Unpolarized states, SU(2) invariant states, and thermal 
states 

The unpolarized states lISTl l32l have isotropic central mo- 
ments for all orders. Due to symmetry, the odd order cen- 
tral moments are identically zero. For an A^-photon unpolar- 
ized state, the second-order central-moment in any direction is 
(A^) = (§1) ^N{N + 2)/3, and the fourth-order central mo- 
ment is (A4) ^N{N + 2){3N^+6N-4-)/15. 

SU(2) invariant states ll35l constitute a subclass of the un- 
polarized states fn]. They take the form 

E (21) 

where pN is a probability distribution and Ia' is the iden- 
tity operator in the Mh manifold. A simple example of 
a state that is unpolarized but not SU(2) invariant is given 
in Ref . Ii32l . In Table |IV] we have tabulated the measured 
data for the unpolarized (and SU(2) invariant) two-photon 




(b) /T ■ ■ 




FIG. 3: The second-order polarization central moment (A„) for the 
state (|2,0)(2,0| + |0,2)(0,2|)/2. Theoretical plot in (a) and the ex- 
perimental results in (b). 

state (|2,0)(2,0| + |l,l)(l,l| + |0,2)(0,2|)/3. The measured 
Stokes vector for this state is (-0.07,-0.10,0.01). As can 
be read essentially directly from the table, the experimentally 
obtained function (A^) is more or less a sphere. The exper- 
imentally obtained extremal values of the function are 2.72, 
2.68, and 2.62, indicating an error due to fluctuations of about 
±2%. 

The thermal states, finally, constitute a subclass of the 
SU(2) invariant states with = iV"/(l +N)^+\ where 
is the average excitation [e\p{hv /kT) — 1]^', where h is 
Planck's constant, v is the optical frequency, k Boltzmann's 
constant, and T the temperature. Their second- and fourth- 
order central moments are given above, and their excitation 
probability averaged second- and fourth-order central mo- 
ments are + 2N^/3 and + 26A^V3 + 12A^^ + 2W'*/5, re- 
spectively. 



E. A two-photon mixed state 

By mixing states, one can get quite complicated pola- 
rization characteristics. Below, we shall elaborate on this 
for three photon states. For two-photon states the parame- 
ter space is of course smaller. For example, for the state 
(|2,0) (2,0| + |0,2) (0,2|)/2, the second-order central-moment 
(A^)2 depicted in Fig. bI (theory on the left, experiments 
on the right). The measured Stokes vector of this state is 
(—0.11, —0.10, 0.00) . To better appreciate the random and the 
systematic errors of the measurement, a cut through Fig.[3]in 
the (A^) = plane is shown in Fig. U\ Again it is seen that 
the experimentally obtained figure is not aligned with the in- 
tended orientation, but that the experimentally obtained figure 
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FIG. 4: (A^) in the plane defined by (A^) = for the state 
(|2,0)(2,0| + |0,2)(0,2|)/2. The theoretically expected resuU is 
drawn solid (blue), the experimentally obtained result is drawn 
dashed (red), and the experimental result, solidly rotated so that its 
eigenvectors coincide with the intended eigenvectors, is drawn dotted 
(black). 



is rotated about 10 degrees around the Si axis. The figure's 
extremal values are 3.99, 2.03, and 1.97. Reorienting the fig- 
ure by a solid rotation so that the measured and intended axes 
coincide, one obtains the dotted curve in Fig. |4] By so re- 
moving the systematic errors, one obtains a very good (±2%) 
agreement between the experiments and the theory. 



F. Some 3-photon states 

Here we explore the polarization characteristics of a few 
states up to the third order and show that for N = 3 there ex- 
ist six classes of states if they are sorted according to their 
first, second, and third-order polarization central moments. 
Before giving examples of the classes, it is helpful to retain 
the uncertainty relation Q. In the third excitation manifold 
this means that the sum of the second-order polarization vari- 
ances must lie between the values 6 and 15. In order to have 
an isotropic second-order central moment, the diagonalized 
covariance matrix F should be proportional to the 3x3 unit 
matrix, and the relation above dictates that the proportional- 
ity factor must be in the range between 2 to 5. In fact, only 
minimum sum uncertainty states will reach the lower limit in 
Q and such states have an anisotropic second-order polariza- 
tion central moment. We conjecture that the lower limit for 
an isotropic second-order central moment is in fact (So) so 
that the minimum uncertainty sum for second order isotropic 
states is 3{Sq) (and specifically 9 for three-photon states). The 
corresponding state is |A^,0)(A^,0|(1 ± [{N - 1}/nY^^)/2 + 
|0,A^)(0,A^|(1 T [{A^- l}/A^]'/^)/2, and in this specific mani- 
fold (l/2±6-V2)|3^0)(3,0| + (1/2t6-i/2)|0,3)(0,3|. 

Since states that lack first-order polarization but are sec- 
ond order polarized have already been discussed, we shall 



now look at states that have an isotropic, second-order po- 
larization central moment, but that may have higher-order 
polarization-structure. Applying the requirements for a state 
to have isotropic polarization up to second order, one can de- 
rive such a three-photon state's density matrix p to be of the 
form 



P = 



/ Pii 

P*n 
P*n 
\P*u 



Pl2 

l-3pii 



Pl3 

-V3pi2 
3pii-5 

Pl*2 



Pl4 
-Pl3 
Pl2 
-Pll / 



(22) 



Here, pi i is real and subject to the restriction l/6<pii<l/3, 
whereas pi2, pn, and pi4 may be complex. Of course, the 
general coherence property for the off-diagonal coefficients 
\Pjk\ < \/PjjPkk holds and imposes additional (but simple) re- 
strictions on the matrix, once one has chosen pn. 

Note that the matrix above defines the sufficient conditions 
for a density matrix to have vanishing first and second-order 
central moments, but that it does not include all necessary con- 
ditions to make it a density matrix. That is, it may be that, 
for certain choices of parameters, the matrix is not strictly 
non-negative. Hence, the reader is warned that when using 



Eq. ( 22 1, to make sure that the ensuing matrix is non-negative 



For p of the form in Eq. (22 1 one finds that (8)3 = (0,0,0), 
that (A?) 3 = 5 for i ~ 1,2,3, and hence that 



(23) 




One can also deduce from ( 22 1 that there is no pure, three 



photon state that is unpolarized to second order. This follows 
from the condition that for a pure state, \Pjk\^ = PjjPkk- Ap- 



plied to (22 1 one gets the three conditions |pi2| — Pii(l ^ 
3pn), 3|pi2p = (1 - 3pii)(3pii - i), and \pn\^ - (3pn - 
5) (5 ^Pii)- The first two of these equations demand that 
Pi 1 = 1/3 — > P12 = 0, but this value does not satisfy the third 
equation. 

An already discussed class of states are the unpolarized 
states. This is the smallest class of three -photon states, be- 
cause there is only one such 3-photon state. The state has, of 
course, isotropic polarization properties of all orders, but re- 
quiring this property for only the lowest three orders uniquely 
singles out this state. 

The mixed state ||3,0)(3,0| + i|l,2)(l,2| + i|0,3)(0,3| 
lacks first-order polarization, has (A^)3 = 5, but has third- 
order polarization structure. Its third-order polarization cen- 
tral moment is shown in Fig. [5] This state is thus unpolarized 
to second order. 

The mixed state (|3,0)(3,0| + |0,3)(0,3|)/2 has vanishing 
first and third-order central moments in all directions, but has 
an anisotropic second order polarization central moment, with 
the predominant fluctuations along the ^3 axis. (A^)3 has a 
"peanut" shape (similar to Fig. [3]l with "semi-axes" lengths 
(Aj)3 = (A2)3 = 3 and (A3) 3 = 9. This is thus a maximum 
uncertainty state. 

The pure state (|0,3) + |3,0))/V2 lacks first-order polari- 
zation, has (A2)3 identical to the (|3,0) (3,0| + |0, 3) (0, 3|)/2 
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FIG. 5: The absolute value of (A^)3 for the state =^|3,0)(3,0| + 
i|l,2)(l,2| + i|0,3)(0,3|. 
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l(Ai>i 2 

FIG. 6: The absolute value of (A^)3 for the state (|0, 3) + 13,0))/\/2. 

mixed state, and its third-order polarization central moment is 
shown in Fig.|6] It is also a maximum uncertainty state. 

Changing the mixing ratios somewhat, one finds that the 
mixed state ^|3,0>(3,0| + i|l,2)(l,2| + ^|0,3>(0,3| also 
has no first-order polarization, but second- and third-order 
structure. The second-order central moment is again "peanut 
shaped" with "half-axes lengths" (19, 13, 13)/3, and (A^)3 is 
similar to that in Fig. |5] Hence, comparing this state with 
the state (|0, 3) + |3,0))/ \/2 they have very similar polariza- 
tion properties in the first two orders (and both are maximum 
uncertainty states), but their third-order properties are vastly 
different. 

Finally, an example of a state that has first and third-order 
polarization structure but has an isotropic second-order cen- 
tral moment is the mixed state 5||3,0)(3,0| + 55|1,2)(1,2| + 
j^|0,3)(0,3|. This state has the Stokes vector (-1,0,0), 
(A2)3 = 14/3, and (A^)3 is shown in Fig.jv] 

It is not possible to find three-photon states that have 
first-order polarization but a vanishing third-order polariza- 




FIG. 7: The absolute value of (A^)3 for the state |||3,0)(3,0| + 
i||l,2)(l,2| + i|o,3)(0,3|. 



tion central moment. The above Ust of possible polarization 
classes exhausts all the combination of polarization structures 
up to third order and shows that six different classes exist 
out of the total of eight a priori possible combinations. The 
classes, with example of associated states, are tabulated in Ta- 
ble H] 



IX. CONCLUSIONS 

We have developed a systematic method, using central mo- 
ments, for assessing the polarization characteristics of quan- 
tized fields. The method goes well beyond the "standard" 
method that only considers first-order moments, and that 
moreover, averages over the excitation manifolds. We have 
shown that there exist a rich "zoo' of polarization states, in- 
cluding, e.g., states that are unpolarized up to a given order 
but that have higher-order structure (so called hidden pola- 
rization). However, as expected, for most states the polari- 
zation characteristics are dominated by the first and second- 
order behavior, as higher-order polarization moments always 
contain "beating terms" originating from lower orders. Some 
states, however, show polarization structure that is dominated 
by higher-order moments, and examples of such states are 
given. 

The suggested method is not the only way to fully charac- 
terize the polarization of quantized fields. In particular two 
more-or-less equivalent methods are mentioned, namely gen- 
eralized coherence matrices 1 19, J3J and the expectation val- 
ues of all combination of Stokes operators ||32]| . 

Since a state is not fully specified by its polarization prop- 
erties, it comes as no surprise that polarization tomography is 
less resource demanding than full state tomography. We have 
quantified this difference and indicated a "recipe" for deter- 
mining polarization properties up to a certain order 
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Appendix 

The experiments were performed by using spatially non- 
degenerate, photon-pair states generated in the process of 
spontaneous parametric down-conversion. The photon pair 
centered at 390 nm was generated in a 2 mm thick type-1 )3- 
barium-borate (BBO) crystal pumped by a femtosecond laser 
pulse centered at 780 nm wavelength. The photon-pairs were 
subsequently filtered by an interference filter with a 4 nm full 
width at half maximum (FWHM) bandwidth. The photon 
pair were brought to the inputs of Hong-Ou-Mandel (HOM) 
interferometer |36:|. When the photons' wavefunction over- 
lap in the HOM interferometer, either the state or the 
state |2,0) can be postselected, dependent on the relative po- 
larizations of the incident photons. The generation setup is 
described in more detail in |,37J . 

In order to measure the first and second order Stokes pa- 
rameters, a polarizing beam splitter (PBS) is positioned at the 
measurement stage. The measurement basis is changed by 
means of one half- and one quarter- wave plate which are set in 
front of the PBS. At each output of the PBS, a two-photon de- 
tector is simulated by a 50:50 fiber beam splitter (FB) and two 
single-photon detectors (PerkinElmer, SPCM-AQRH). The 
relative coincidence detection-efficiencies are estimated from 
the FB transmittance. Subsequently, the photon detection- 
efficiency of each single-photon detector-channel is used to 
calibrate the measurement of the Stokes parameters. The 
relative coincidence detection efficiencies of the four detec- 
tors are 0.91 : 0.91 : 0.82 : 1, for |1, 1) and 0.75 : 0.76 : 1 : 
0.57, for |2,0). In order to achieve full information about the 



first and second order Stokes parameters, we measured these 
coincidences in six distinct measurement bases. For a precise 
measurement, we measured coincidences three times and each 
measurement is done for 3 s. The central moments are ob- 
tained from the measured Stokes operator at the six different 
directions, and then solving the equation system generated by 



Eq. (111. In the tables, averages of the measurement and the 
estimated errors due to fluctuations are presented. As can be 
seen, what looks like errors far exceeding the error bars (e.g., 
for (52^3 +53^2) of state |2,0)) are actually systematic errors 
due to imperfect polarization optics. These can be corrected 
by properly aligning the experimental Poincare axes with the 
"theoretical" axes. When this is done, as is shown in Figs. [2] 
and|4j the experimental figures are coinciding to within a few 
percent with the theoretical ones, showing that the errors due 
to fluctuations are relatively modest and within the relative 
range indicated by the error Umits in the tables. 
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Polarization transformation invariant 


State 


{Sn)N 






i/4 


Yes 


Yes 


Yes 


i|3,0)(3,0| + i|l,2)(l,2| + i|0,3)(0,3| 


Yes 


Yes 


No 


2(|3,0)(3,0| + |0,3)(0,3|) 


Yes 


No 


Yes 


(|3,0) + |0,3))/x/2 


Yes 


No 


No 


||3,0)(3,0| + i5|i,2)(l,2| + T^|0,3)(0,3| 


No 


Yes 


No 


|3,0) 


No 


No 


No 



TABLE I: A table of states exemplifying the six different 3-photon polarization classes. 





Theory 


Experiment 


Operator 


1 St order 


2nd order 


1 St order 


2nd order 


(Si> 





2 


-0.19±0.06 


2.06 ±0.03 


(Si) 





2 


0.12±0.04 


2.04 ±0.02 


{S3) 


2 


4 


1.97±0.01 


3.93 ±0.02 


{S1S2+S2SI) 









-0.09±0.10 











0.54±0.08 











-0.13±0.15 



TABLE II: The experimental data for the state |2,0). 





Theory 


Experiment 


Operator 


1st order 


2nd order 


1st order 


2nd order 


(si> 





4 


-0.01 ±0.04 


3.98 ±0.00 


(S2> 





4 


-0.08 ±0.03 


3.93 ±0.03 










0.01 ±0.02 


0.15 ±0.05 


{SiSi+SiSi) 









0.03 ±0.02 











-1. 67 ±0.14 


{SjSi+SiS}) 









0.12±0.16 



TABLE ni: The experimental data for the state 1 1, 1). 





Theory 


Experiment 


Operator 


1 St order 


2nd order 


1 St order 


2nd order 


(si> 





8/3 


-0.07 ±0.03 


2.69 ±0.03 







8/3 


-0.10±0.02 


2.68 ±0.03 


(S3> 





8/3 


0.01 ±0.01 


2.67 ±0.02 


(S1S2 + S2S1) 









0.00 ±0.07 


{S2S3+SA) 









-0.09 ±0.06 


(S3S1 +SiS}) 









0.04 ±0.09 



TABLE IV: The experimental data for the state (|2,0)(2,0| + 
|l,l)(l,l| + |0,2)(0,2|)/3. 





Theory 


Experiment 


Operator 


1 St oi'der 


2iid order 


1 St order 


2ird order 


..•''Y; 







0,1 1 ^0,04 


2.04:1:0.04 


{Sl) 





2 


-0.10 ±0.03 


2.05 ±0.04 







4 


0.00 ±0.01 


3.93 ±0.02 


{S1S2+S2S1) 









-0.01 ±0.11 


{S2S3+S3S2) 









0.69 ±0.07 


(SiSt +S1S3) 









-0.01±0.11 



TABLE V: The experimental data for the state (|2,0)(2,0| + 
|0,2)(0,2|)/2. 



